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Abstract 

This paper, about a fluid-like system of spatially confined particles, reveals the analytic 
structure for both, the canonical and grand canonical partition functions. The studied 
system is inhomogeneously distributed in a region whose boundary is made by planar faces 
without any particular symmetry. This type of geometrical body in the <i-dimensional space 
is a polytope. The presented result in the case of d = 3 gives the conditions under which 
the partition function is a polynomial in the volume, surface area, and edges length of the 
confinement vessel. Equivalent results for the cases d— 1,2 are also obtained. Expressions 
for the coefficients of each monomial are explicitly given using the cluster integral theory. 
Furthermore, the consequences of the polynomial shape of the partition function on the 
thermodynamic properties of the system, away from the so-called thermodynamic limit, 
is studied. Some results are generalized to the <i-dimensional case. The theoretical tools 
utilized to analyze the structure of the partition functions are largely based on integral 
geometry. 
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I. INTRODUCTION 



Thermodynamic properties of fluids are relevant to biology, chemistry, physics and en- 
gineering. In most cases of interest the fluid system is inhomogeneously distributed in the 
space, confined to a region of finite size and constituted by a bounded number of molecules. 
A typical example of this kind of systems is that of fluids confined in pores. They consti- 
tute a prototypical inhomogeneous system which occupies a small region of the space and 
may involve a small number of particles. It is well known that the spatial distribution of 
a fluid confined in a small pore may follow or not the symmetry of the cavity , and 
that its properties may be strongly influenced by the geometry of the container. To study 
the thermodynamic properties of these systems it is customary to introduce several strong 
assumptions or approximations that simplify the analysis. An usual approach is to treat 
the system as it were homogeneous (which necessarily implies that the system completely 
fills the space, and thus, involves infinitely many particles). A second usual approach as- 
sumes that the inhomogeneous fluid is spatially distributed in several regions, each one with 
homogeneous properties, while the inhomogeneous nature of the system is concentrated in 
regions with vanishing size (surfaces, lines and points). Frequently, this scenario is comple- 
mented with the assumption that the spatial distribution of the fluid involves continuous 
translational and/or rotational symmetries. Finally, a third approach only assumes that 
the inhomogeneous fluid takes spatial configurations with those symmetries. In general, 
the assumed symmetric distribution of the fluid could be attained spontaneously (as in the 
case of free drops or bubbles) or could be induced by an external potential as it is the case 
of confined system that are constrained to regions with simple symmetry: a semi-space, a 
slit, an infinite cylinder, or a sphere (for example in wetting and capillary condensation, 
phenomena). Anyway, the continuous symmetries (translations and/or rotations) play an 
important role providing the bases to identify the extensive and intensive magnitudes which 
enable the development of thermodynamic theory [3|. 

Moreover, thermodynamics only provides an incomplete set of relations between intensive 
and extensive magnitudes which must be complemented with other sources of information to 
obtain the thermodynamic properties of a given fluid system. In any case, these sources, that 
may be experimental, theoretical, or based on numerical simulation, also involve assump- 
tions or approximations related to the existence of continuous symmetries. In this sense, 
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our statistical mechanical and thermodynamical approach to the study of inhomogeneous 
fluids appears to be intrinsically entangled with some hypothesis about the constitution and 



behavior of the system under study 



2J|. Particularly, I refer to three hypothesis concerning 



the fluid system: the symmetry of its spatial distribution, the large volume occupied and the 
large (and unbounded) number of particles involved. Of course, one may ask about if those 
hypothesis are or not a central part of the theories. What happens with inhomogeneous 
fluid-like systems which do not necessary attain spatial distribution with simple symmetry 
and/or occupy small regions of the space and/or are constituted by a small number of par- 
ticles? Can we apply statistical mechanics and thermodynamics to study their equilibrium 
properties? How can we do that? The analysis presented below attempts to advance in 
the understanding of these questions. From a complementary point of view, this work also 
deals with the long standing objective of finding fluid-like systems that are exactly solvable, 
i.e., its partition function integral can be integrated and thus transformed into an analytic 
expression. Each of these systems provides the unique opportunity of testing some of the 
fundamental hypothesis of the theoretical framework. 

This paper is devoted to analyze, using an exact framework, few- and many-body fluid 
like systems confined in cavities without continuous symmetry. In particular, we focus on 
their partition function and thermodynamic magnitudes as functions of the spatial set that 
defines the region in which particles are allowed to move. In fact, it is shown that under 
certain conditions the partition function of the fluid-like system is a polynomial function 
of certain geometric measures of the confining cavity, like its volume and surface area. 
The adopted approach is theoretical and exact, being our core results largely based on 
integral geometry. The rest of the manuscript is organized as follows: In Section |TT] it 
is presented the general statistical mechanics approach to canonical and grand canonical 
ensemble with two modifications: the cutoff in the maximum number of particles which 
enables the analysis of few-body open systems and the generalization of cluster integrals 
to inhomogeneous fluids. Sections IIHI and IIVI are devoted to study the structure of cluster 
integrals with particular emphasis in the case of a polytope-type confinement. There, two 
theorems and their corollaries are demonstrated, which constitute the main result of the 
present work (PW). Extensions to other type of confinements are discussed in Section IVl 
while the consequences of the cluster integral behavior in the thermodynamical properties 
of the system are studied in Sec. IVIl The final discussion is given in Section IVHl 
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II. PARTITION FUNCTION OF AN INHOMOGENEOUS SYSTEM 



In this work we consider an open system of at most M particles that evolves at constant 
temperature in a restricted region A of the d- dimensional euclidean space M, d . The region 
A is the set of points where the center of each particle is free to move and has a boundary 
dA. From here on the system will be shortly referred to as the fluid while A will be referred 
to as the region or set where the fluid is confined. In fact, the system is an inhomogeneous 
fluid due to the existence of dA. The restricted grand canonical partition function of the 
fluid is 

M 

Sm = $>^-, (1) 

j=0 

where M is the maximum number of particles that are accepted in A. The absolute activity 
is x = exp(/3/i), being fi the chemical potential while /3 = 1/ksT, ks and T are the 
inverse temperature, Boltzmann constant and temperature, respectively. Em is a M-degree 
polynomial in x and its j-th coefficient is the canonical partition function of the closed 
system with exactly j particles confined in A. The grand canonical partition function of the 
fluid without the cutoff in the maximum number of particles can be obtained from lim Sm- 

M -toe 

On the other hand, the canonical partition function of the system with j particles is 

Qj = IjA-^Z, , (2) 

where Ij = is the indistinguishably factor, A = h/(2irm hsT) 1 / 2 is the thermal de 
Broglie wavelength, m is the mass of each particle and h is the Planck's constant. Therefore, 
we can transform Eq. ([I]) into 

M 

Z M = J2l j z 3 Z j , (3) 

3=0 

where the activity is z = exp (/3fi) /A d . Finally, the configuration integral (CI) of a system 
with j particles is 

j 

Y[e k Yl ei n d j r, (4) 

k=l <l,n> 

being Z = 1 [alternative dimensionless definitions for z and Zj may be obtained by in- 
troducing the volume of the system or the characteristic volume of a particle in Eq. (j4j)]. 
Here, e; n = exp [— /3(f)(ri n )] is the Boltzmann factor related to the spherically symmetric pair 
potential <fi between the I and n particles, which are separated by a distance ri n . It will be 
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assumed that has a finite range £, being 0(r; n ) = and e/ n = 1 if r/ n > £. This assump- 
tion is not very restrictive because any pair interaction potential can be approximated by 
truncation at a finite range, e.g. it is frequent to study the Lennard Jones fluid cut and 
shifted at r = 2.5a % |5fl- The indicator function 

{1 if r k e A , 
( 5 ) 
if r k £ A , 

is the Boltzmann factor corresponding to the external potential produced by a hard wall 
confinement. Note that, the integration domain in Eq. (j3j) is the complete space due to the 
spatial confinement of the particles in A is considered through ej.- The CI is a function of 
the set A and a functional of <p. Z{ itself is given in its full generality by [6fl 



IT— C? 



m,l \i\ 

i=l 1 



(6) 



where the sum is over all sets of positive integers or zero m = {mi, 777.2, • • •} such that 
YjI=i = 3- Eq- © shows that Zj is a polynomial in Ti, ... ,Tj which are essentially the 
(reducible) Mayer cluster integrals for inhomogeneous systems. This equation is obtained 
through the following procedure: replace in Eq. (j3J) each e;„ using the identity e\ n = 1 + f\ n 
(this Eq. defines //„,), distribute the products, and collect all terms of the integrand which 
consist of groups of particles that conform a cluster, in the sense that they are at least 
simply connected between them by /-functions. Note that our assumption about e\ n = 1 
for ri n > £ implies f] n = if r^ n > £ showing the spatial meaning of the cluster term. The 
procedure gives Eq. Q with 

n= [■■■[ S h2 _ tl f[e k d l r , (7) 

k=l 

Si,2,... ti = J2 Uf - ( § ) 

cluster <l,n> 

where Si 2,. * is a sum of products of fi n functions with 1 < I, n < i that involves all the 
products of / functions, which can be represented as a connected diagrams (clusters) with i 
nodes and fi n bonds. Clearly, 73 depends on T, and A. If one assume that A is very large 
the usual homogeneous system approximation gives n — > HVbi, where V is the volume of A 
and 6jS are the Mayer cluster integrals which depend on T and (p. 
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The inversion of Eq. ([6]) gives the following expression for the dependence of Tj with the 
CIs (Eq. (23.44) in Ref. |g) 



n 



1 (Zj 



Ujl \j\ 



(9) 



where the first sum is over all sets n = {rii, ri2, ■ • •} with nj non- negative integers such that 
Yjj=i3 n j = i an( ^ r i = ^i^i w ith &i = 1. From Eq. Q it is apparent that Tj depends on 
Zj with j = 1, . . . , i, and thus, the TjS with i — 1, ...,k and the Z^s with j — 1, . . . , k involve 
the same physical information. On the other hand, one can return to the Eqs. (JT]) - (j3)) to 
observe that they can be re-written in this alternative form: replace M — > oo in Eq. ([1]) but 
assume ZM+k = with k — 1, . . . , oo. In this context Eq. (Q shows that cluster integrals Tj 
with i = 1, M are not affected by the restriction j < M, being Tm+i the first affected r 
(because it depends on Zm+i, which is zero). 

Before ending this section we wish to focus on a relevant characteristic of H^f, Zj and Tj 
functions. Let us define § = |.A/«4 C M d j (the set of all the subsets of M. d ), for fixed T, 
and z one can write Sm(A) : § — > K which implies that Ha/(-4) may depend on the shape 
of A. Clearly, the same argument applies to Zj and Tj which may also depend on the shape 
of A. Here we anticipate the principal result of PW, related with this non-trivial shape 
dependence, that will be demonstrated in the Sees. IIIII and HVl Thus, we turn the attention 
to a fluid confined by a polytope A. For simplicity we focus in the three dimensional case, i.e. 
a fluid confined by a polyhedron. If a system of particles that interact via a pair potential 
of finite range £ is confined in a polyhedron A such that its characteristic length £(A) [see 
Eq. (125]) ] is greater than kC, + C (being C a constant) for some integer k > 1. Then the i-th 
cluster integral Tj with 1 < i < k is a linear function in the variables V, A and {L\, L2, ...} 
(the length of the edges of ^4). In fact 

Ti /i\ = V b l -Aa l +Y, LncU + £ c\ n , (10) 

nedges nvertex 

where the coefficients bi and <2j are independent of the shape of A while c e in and cj n are 
functions of the dihedral angles involved. Besides, all the coefficients depend on the pair 
interaction potential and temperature. Expressions similar to Eq. fllOp are also found for 
the euclidean space with dimension 2 and 1, while they are conjectured for dimension larger 
than 3. Two non-trivial consequences derive from the Eq. ( ]T0]) . On one hand it implies 
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that, if the system of particles confined by A involves N particles and £(A) > iV£ + C, then 
Zn is polynomial on V, A and {L\, L 2 , ...}. On the other hand it implies that, if the system 
of particles confined by A is open, involves at most M particles and £(«4) > M£ + C, then 
Sjv/ is a polynomial in z, V, A and {L 1; L 2 , ...}. 

III. THE PROPERTIES OF SOME FUNCTIONS RELATED TO 77 

In this section we analyze the properties of some many-body functions related to the 
partial integration of the cluster integral r< [Eqs. (j7|) and (jSJ)]. This analysis will be comple- 
mented in the next section where we will reveal the linear behavior of n. In the following 
paragraphs several definitions are introduced and two different proofs of the locality and 
rigid invariance of functions related to the partial integration of n are presented. Both 
proofs are necessary for clarity. In the first approach the mentioned properties of those 
functions are demonstrated and it is found a cutoff for its finite range, while in the second 
approach (which is more complex than the first one) the properties are proved and a better 
bound of the finite range is obtained. 

We define, G(A, r) : § x M. d — > R is a local function in r over A with range A (from 
hereon local with range A) if its value is entirely determined by the set U (r, A) D A that is 
G(A, r) = G(U (r, A)n^4) where U (r, A) is the ball centered at r with radius A. This definition 
can be generalized to functions of several variables in the following way: G(A, ri, . . . , r n ) : 
§ x M. d * n — y R is said to be a local function in over A with range A if G(A, ri, . . . , r n ) = 
G(U (ri, A) n A, ri, . . . , Ti-i, r i+x , . . . , r n ) for {ri, . . . , rj_i, r i+ i, . . . , r n } C U (r h A) n A and 
G(A, Ti, . . . , r n ) = when exist U (r^, A) fl A. 

Let G(A, r) be a local function with range A. We say that G(A, r) is invariant under rigid 
transformations (the elements of the euclidean group, i.e., any composition of translations, 
rotations, inversions and reflections) if V rigid transformation R, R[U (r, A)fb4] = U (r', X)C\A' 
with r' = Rr implies G[U(r, A) DA] = G[U(r', A) PI A']. The generalization to functions with 
many variables is: let G(A, rj, . . . , r n ) : § x IR d * n — y R be a local function in with range A, 
we say that G(A, T\, . . . , r n ) is invariant under rigid transformations if V rigid transformation 
R, R[U(r i} X)nA] = Ufa, A) n A' implies G[U(r u A) n A, v x , ..]= Gp^, A) n A, r' 1: . . .], 
with r'j = Rrj. A direct consequence of the locality with range A and rigid transformation 
invariance of a bounded function G(A, r) is that it attains a constant value for all r G A 
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such that U(r, A) n OA. 

Theorem 1: Let ejt(r) and Si j2l ...,i be the Boltzmann factor and cluster integrand intro- 
duced in Eqs. (jSJ) to (jSJ), and A € § a set in IR d for which the integral in Eq. (II ip is finite. 
Then 

E 1 (A,r 1 )= l~( 5 ljW J^e fc dr 2 ---dri (11) 

J J k=2 

is a local function with finite range invariant under rigid transformations. 
First Proof: Consider the following 2-body function 

Ei(r u y 2 , ...,y») = 5*1,2,.. „i , (12) 

with rj the position of an arbitrarily chosen particle, and the rule to obtain the (n — l)-body 
function from the n-body one given by 

S n _ 1 (^4,r 1 ,y 2 ,...,y„_i) = J E n (A,r 1 ,y 2 , . . . ,y n ) e n dy n , (13) 

where 2 < n < i and yj = Xj — r x is the coordinate of particle j with respect to particle 1. 

Si,2,...,i has range 

(14) 

being Si^,. i = if r a b = \r a — r^l > (i — 1) £ for at least one pair of particles a and b 
in the cluster. This property derives from the fact that Si 2,...,t contains the open simple- 
chain cluster term S' = /12/23 • • • fi-i,i that enables that two particles reach the maximum 
possible separation (z — 1) £ of all the cluster terms in Sip i (there are i\/2 terms of this 
type). Even more, one can show that S\z i is local with range Naturally, Eq. (fT2l) 
shows that Ei(ri, y 2 , . . . , y«) is also local with range q. The integration in Eq. ( fT3l) applied 
to Ei(ri, y 2 , . . . , yi) implies that Ei_i(A, r 1; y 2 , . . . , yj_i) is local in ri with range q. We 
proceed by induction. Let us assume that for some n < i, E n (A, Ti, y 2 , . . . , y n ) is local 
with range <j, then by Eq. (TT3]) E n -i(A, T\, y 2 , . . . , y n -i) is also local with the same range. 
The procedure continue until the function Ei(A, 1*1) is reached. Therefore, we find that the 
E n (A, Ti, y 2 , . . . , y n ) functions with 1 < n < i are local in ri with range <; [even more, we 
have obtained that E n (A, ri, . . . , r n ) are local in rj with range for any 1 < j < n\. 

Taking into account that i?j(ri,y 2 , . . . , y$) has range q and that e, = ej(*4., r-j) [see Eq. 
(jSJ)] the right hand side of Eq. 0131) for n = i can be written as 

E i (ri,y 2 ,...,y i )dy i . (15) 



I 
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It is convenient to express the coordinates (ri, y 2 , . . . , y») in terms of the rigid transformed 
coordinates (r'^y^, . . . , yQ (related each other by rj = and y^ = R~ l y'j = R~ lr 'j ~ 

i2 -1 ri with R- 1 the inverse of R) and to change the integration variable to y- (note that 
the Jacobian is one). Thus we find 

f ^(irM, i2-y 2 , . . . , iry ) ^ , (i6) 

7R[C/(r ll? )n^] 

where the integration domain is equal to [/(r^q) fl A' by hypothesis. Given that Ej, is 
invariant under any rigid transformation applied to the coordinate of the particles one can 
drop each R~ l in Eq. 0161) and return to the original form introducing ej = ei(A', r-) 

I E i (r' l ,y' 2 ,...,y' i )e i dy' i , (17) 

which is the definition of _E i _ 1 ( v 4 / , r' 1; y' 2 , . . . , y^j). This shows that E 1 ^^, 1*1, y2> • • • , yi-i) 
is invariant under rigid transformations. Again, we proceed by induction. Let us assume 
that E n (A, y2, . . . , y n ), a local function in r x over ^4. with range q, is invariant under 
rigid transformations. Taking into account that e n = e n (A, r„) for the right hand side of 
Eq. (TT3)) we can write an expression similar to Eq. 015]) but replacing £?j(r 1; y 2 , . . . , y*) by 
E n (U(ri, <;) fl ^4, y2, . . . , y n ) and % by n. Turning to transformed coordinates we find 

f EniRT^Uft, n .A'], i2-y 2 , • • • , i2-y n ) dy' n , (18) 

>/R[[/(ri,?)n.A] 

where the integration domain is equal to U(r[,q) fl A' and we used that U(ri,q) PI A — 
-R _1 [t/(r' 1 , <j) fl*4/]. Given that E n is invariant under rigid transformations we can drop each 
R^ 1 in the arguments of E n in Eq. ([TBI , use the finite range of E n to split the argument 
U(r[,q) fl A' into (A.', r^) and introduce e n (A',r' n ) to obtain 



/ 



E n (A',r , v y 2 ,...,y , n )e n dy , n , (19) 



which is the definition of E n _ x {A', r' 1; y' 2 , . . . , y^_ x ). Therefore, £ n _! (A, r 1; y 2 , . . . , y„_i) 
is local of range q and rigid transformation invariant. The procedure continue until the 
function E\ is reached. Therefore, returning to the original coordinates, we obtain that 
E n (A, Ti, . . . , r n ) with 1 < n < i is local in any Tj (1 < j < n) with range q and rigid 
transformation invariant. In particular, for n — 1 it implies that Ei(A, ri) is local in r x with 
range q and rigid transformation invariant. ■ 
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Second Proof: This second proof is based on a more subtle consideration about the 
particle labeled as ri in the first proof. We introduce the coordinate of the central particle 
of the cluster, r c , and the relative coordinates yj = rj — r c of the particle j with respect to r c . 
Again, consider the cluster integral Tj and its simple open-chain cluster S' term. For S' we 
take r c = with I = Integer Part [i/2 + 1] (if i is an odd number r c = r( i+1 )/ 2 is the position 
of the middle-chain particle while if i is an even number r c = rj/2+1 is the position of one 
of the pair of particles that are at the middle of the chain). We note that S' 12 j is zero if 
|r a — r c | > (/ — 1) £ for at least one particle a in the cluster. The procedure to obtain r c for 
all the other terms in S\ : 2,...,i is as follows: for a given cluster 5"' of i particles take iteratively 
each pair of particles, separate them to find the maximum possible elongation distance 
under the condition S" 7^ 0. Let the more stretchable chain of particles (that could be 
non-unique) be a chain of k particles with end-particles a and b. Thus, |r — r&| < (k — 1) £ 
with 2 < k < i and for this cluster we can define r c = with I = IntegerPart [k/2 + 1]. By 
using this second approach we find that Si^,...^ has finite range 

^ = ({ - 1) £/2 if % is odd, q = i£/2 if i is even . (20) 

Once r c is identified for each cluster term in S\ t 2,...,i we can rename r c as ri and follow the 
procedure developed in the first proof of the theorem. ■ 
Based on Eq. ( 120]) along with the assumption that the range of E(r) must be a unique 
function of the maximum elongation length of Si,2,...,i independently of the parity of i we 
have the following guess 

q = (i- i)f/2 Vi e N. (21) 

We may mention that Eq. (12T1) can be demonstrated for the case of A being a convex body 
[by virtue of Eq. (T2"0"j) one must focus on the case of even i\. The demonstration follows a 
procedure similar to that used to obtain the Eq. (l20l) . but for the case of even i one define 
r c = (r-j/2 + r(i/2+i)) /2. The development made in this section enable to write Eq. (J7|) as 

t 1 = J E{x x )e x dxx , (22) 

with two different definitions for E{ji). Both definitions and any other possible approach 
must give mathematically equivalent expressions for Tj. 
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IV. INTEGRATION OVER A POLYTOPE-SHAPED DOMAIN 



In order to analyze the implications of the domain's shape on certain type of integrals 
it is necessary to introduce some notions about sets B, C C R d . They are summarized in 
this and the next paragraphs. The closure of B is cl (B) = B U dB. It is said that B is a 
closed set if B = cl (B). Besides, the interior of B is int (B) = B \ dB. It is said that B is an 
open set if B = int (B). There are sets that are neither closed nor open. {B±, . . . , B n } is a 
partition of the non-empty set B if B = Ufe=i &k an d, f° r an h 3 with z 7^ j, Bi 7^ 

and ^ n = 0. A set in R d is connected if it cannot be partitioned in two non-empty 
sets B and C such that cl (B) fl C = i3 fl cl (C) = 0, otherwise it is disconnected. We also 
introduce the concept of connectedness-based partition of a set; a partition {£>i, . . . , B n } of 
B is the connectedness-based partition of B if, for all i, j G N< n with z 7^ j, £>j is connected 
and Bi U i3j is disconnected. The following notions of distances are adopted: the distance 
between the point r G R d and the set B is d (r, B) = min (|r — b| , b G B) with |r — b| the 
usual Euclidean distance between points, while the distance between the sets B and C is 
d (B, C) = min (|b — c| , b G B and c G C). Finally, a set B is said to be convex if every pair 
of points x and y in B are the endpoints of a line segment lying inside B. 

An affine fc-subspace is a linear variety of rank k in R d , whether it contains the origin or 
not. For example, the (d — l)-dimensional affine subspace is a hyperplane. The affine hull 
of a set B C R . , aff is the affine subspace with the smallest rank in which every point 
of B is contained. Given that d is a fixed parameter from hereon we will refer to the rank 
of aff (B) as the dimension of the affine space of B, denoted dim [aff (B)} or simply dim (B). 
For B C R d , the sets relcl(S), relint (B) and veld (B) are the relative closure, the relative 
interior and the relative boundary of B, i.e., the closure, interior and boundary of B within 
its affine hull, respectively 

Even though a polytope is essentially a mathematical entity, in the current work it also 
has a physical meaning because it serves to characterize the shape of the vessel where 
particles are confined. This duality makes difficult to start with a simple and satisfactory 
geometrical description of the kind of polytopes that are relevant for our purposes. We 
begin a somewhat indirect approach that ends in the definition of the set of polytopes that 
have simple boundary (a well behaved one), in the sense that the boundary can be dissected 
in its elements or faces. To make further progress it is convenient to introduce the convex 
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polytopes which, from a geometrical point of view, are simpler than the general polytopes. 
A convex polytope in M. d is any set with non-null volume given by the intersection of finitely 
many half-spaces Q (this definition includes unbounded, closed and non-closed polytopes). 
On the other hand, a connected and closed (CC) polytope A is the connected union of 
finitely many closed convex polytopes. We define, the CC polytope A C IR d is a simple- 
boundary (SB) polytope if for each r G OA 3 e G M>o such that VA < A < e, both sets 
int [U (r, A) D A] and int [U (r, A) fl A c ] are topologically equivalent (homeomorphic) to an 
open ball. It is clear that a closed convex polytope is also a CC SB polytope. In this sense, 
we say that the boundary of a CC SB polytope and the boundary of a convex polytope are 
locally equivalent. The SB condition excludes some degenerate or pathological cases e.g. 
a polytope with two vertex or two edges, in contact. On the other hand, the definition of 
SB polytopes includes bounded and unbounded polytopes, non-convex polytopes, polytopes 
with holes or cavities and many kind of faceted knots embedded in M. d . Let F d be the class 
of CC SB polytopes in R d . 

For a given polytope we focus on the partition of its boundary based on its faces. As 
before, we treat first the convex polytope case. Let B G P d be a convex polytope and H a 
hyperplane. If H n dB ^ 0, H n int (B) = and k = dim (H n OB) with < k < d - 1, 
then we say that HndB is a closed fc-face of B [8]. One can demonstrate that the set whose 
elements are all the closed 0-faces and relopen fc-faces with < k < d — l,ofi3isa partition 
of dB. Other polytopes, non-necessarily convex, may admit a similar face-decomposition 
of its boundary. Given a polytope A with face-decomposable boundary, we introduce the 
notation dA m , n (1 < m < d) for the n-th (d — m) -dimensional relopen face of A (from here 
on an open (d — m)-face), while dAd yU is the n-th closed 0-face. The closed 0-faces are the 
vertex of the polytope. The face-based partition of dA is {dAi t \, ■ ■ ■ , dA m ,n, ■ ■ ■ , dAd, n >, ■ ■ ■} 
with dA m ,i H dA m ',j = for every m,i ^ m',j. Note that dA m ,n is [ct4] mn , i.e., the 
(m, n)-element of dA. Besides, we define 



as the set of points in dA that lies in some of the (d — m)-faces of A, with 
{dAm t i, . . . , dA mtn , . . .} a partition of dA m . Furthermore, 




(23) 



n 




(24) 



fit 
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with m going from 1 to d, where {dAi, . . . ,dAd} is a partition of dA. Given any open 
face, the corresponding closed face is its closure. Furthermore, given any closed fc-face with 
k > 0, its relint is the corresponding open A;- face (an open face relative to its affine hull). 
Some A;-faces are designed by their names: a 3-face is a cell (polyhedron), a 2-face is a facet 
(polygon) and a 1-face is an edge (line). 

The next step is to demonstrate that every CC SB polytope is face-decomposable. To 
find the (d — l)-faces of a CC SB polytope one can exploit that locally they are equal to the 
(d — l)-faces of a convex polytope. Let A C ¥ d , r G dA belongs to an open (d — l)-face of A 
if no matter how small is A with A G M> the set int [U (r, A) fl A] is the interior of a half-ball 
(a ball cut by a hyperplane through its center). In this case, relint [U (r, A) fl dA] C dAi 
belongs to a given open (d — l)-face and we introduce the tangent affine space aff (r, dA) = 
aff [U (r, A) fl dA] of this face. Let r, s G dAi with aff (r, dA) = aff (s, dA) and assume 
that exist a path C (r, s) C dAi that connects r and s such that for every x G C (r, s) 
aff (x, dA) = aff (r, dA), then r, s and C (r, s) are in the same open (d — l)-face of A. 
Furthermore, the connectedness-based partition of dAi is the set of open (d — l)-faces of A, 
i.e. {cMi 5 i, • • • , dA\, n , ' ' ' }• Given that dA = U n cl (dAi, n ), thus every r G dA\dAi belongs 
to the relative boundary of two (d — l)-faces (may be more than two). Let dA\^ and dAij be 
two different open (d — l)-faces of A with cl (dAi,i)C]c\ (dAij) ^ 0, then cl (dAi,i)r\c\ (dAij) 
is a closed m-face of A with m = dim [cl (dAi : i) H cl (dAij)]. By analyzing each pair of 
mutually intersecting closed (d — l)-faces one find each of the remaining dA m , n and dA m 
elements with 1 < m < d. Now, given A G P d we can obtain a non-closed SB polytope 
by removing from A one or several of its faces. However, for PW purposes the polytope 
represents the integration domain of a well behaved function. Given that the value of the 
integral is not modified by removing one or several faces, we will only consider the case of 
closed A. 

Let A G ¥ d , we say that the two elements dA m ,i and dA m ',j of dA are neighbors if 
dAm t i C c\(dA m 'j) or dA m 'j C c\(dA mt i), for m,i ^ m',j. On the other hand, they 
are adjacent if cl (dA m ,i) H cl (dA m ',j) ^ 0- It is simple to verify that if two elements are 
neighbors they are adjacent, and the distance between two adjacent elements is zero. We 
define the characteristic length of the polytope A as 

£ (A) = min {d [dA m ,i, dA m ,J 0} , (25) 
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i.e. the minimum distance between two non-adjacent elements of dA. It is clear that if 
A £ F d , then 3 e £ JR> such that £(A) > e. 

In the following lines we state the general proposition for positive integer values of d. The 
cases of d — 1,2, 3 are demonstrated each one in a separated proof. The remaining cases 
that concern every positive integer value d > 4 are not demonstrated here and thus, the 
proposition is for these cases a conjecture. We note that Theorem 2 and its the conjectured 
generalization to IR d given in the following Eq. (126]) strongly resemble the combination 
of Hadwiger's characterization theorem and The general kinematic formula (see Theorems 
9.1.1 and 10.3.1 in pp. 118 and 153 of [9], respectively) and transpires analogies with several 
results of integral geometry. These points will be discussed at the end of Sec. IVIH 

General conjecture (cases d^ 1): Consider A £ P d and let G(A, r) : F d x R d -> E be a 
well behaved function in r £ cl (A) (with fixed A), and a local function with range q £ M>o 
invariant under rigid transformations with 2q < 2, (A) [for fixed A £ P d function G(r, A) is 
simply G(r) : A R]. Then 

t = / G(r)dr 

J A 

d 

= X d co + X A -,ci + X d - m , n c m ,n, (26) 

m=2 ?i m-elem 

where % d is the (^-dimensional measure (Lebesgue measure) of A (i.e. its volume V), X d -i * s 
the (d — l)-dimensional measure of dA, and \ m „ i s the m-dimensional measure of the n-th 
m-element of dA (in particular x 0n = !)• Besides, Co and ci are constant coefficients which 
are independent of the size and shape of A, while c m ^ n is a function of the angles that define 
the geometry of A near its n-th m-element dAd- m ,n far away from its relative boundary. 

Theorem 2 (case d=3): For the case d = 3 the Eq. ( 1261) is 

t = [ G(r)dr 
J A 

= VC + A Cl + L nC2,n + ^ ' ^ 

n edges n vertex 

where V is the volume of A, A is the surface area of dA, and L n is the length of the n-th 
edge. Besides, Cq and C\ are constant coefficients which are independent of the size and shape 
of A, while C2, n is a function of the dihedral angle in the n-th edge, and C3 >n is a function of 
the set of dihedral angles between the adjacent planes that converge to the n-th vertex. 
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Figure 1: Picture describing the partition for A in its principal and skin parts, A and Ag^, 
respectively. 

Proof (case d=3): Consider the set A G P and introduce a partition of A given by 
{^l ,^l s k} with the principal part A = {r G A/d (r, dA) > q} (the inner parallel body of 
.4) and the boundary or skin part A s k = A \ A D characterized by a thickness q. We have 



Here A Q is the region where G(r) takes a constant value. This partition may be obtained by 
wrapping dA through sliding the center of a ball of radius q on dA. Let w (dA, <r) be the set 
of points that are inside of any of this balls, then A Q = A\w (dA, q). Besides, it may also be 
obtained in terms of a Minkowski sum. Let U be the unit ball centered at the origin, then 
the Minkowski sum A c + qU is the outer parallel body of A c [7\ and A — A \ (A c + qlJ). In 
Fig. [U a picture representing this partition for A is shown. There, one can observe dA on 
continuous line and several balls of radius q (in dotted lines) corresponding to the wrapping 
procedure. For the case of A being the shaded region (both the darker and brighter ones), 
the darker shaded region represents A while the brighter shaded one corresponds to 4 s k- 
On the contrary, if A is the non-shaded or white region then A sk corresponds to the white 
region between dA and the dashed line while the rest of the white region represents A Q . For 
both cases a dashed line separates regions A Q and A sk . The integral in Eq. fl27|) gives 



with g(r) = G(r) — cq and Co = G(r) for all r G A . Naturally, g(r) : A — > K is a bounded 
and local function with range ^ invariant under rigid transformations which implies that 
g(r) = g(U(r, q) fl A) and its value is independent of the position and orientation that the 
set U(r, q) C\ A takes on the space, g(r) only depends on the shape of U(r, q) fl A. Besides 
if r ^ A sk then g(r) = 0. 



-sk • 



(28) 




(29) 
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Note that in the context of d = 3 a planar face is simply a face. Consider the set *4 s k 
partitioned in terms of {Ai, A2, A3}, 

A sk = AiU A 2 U A 3 . (30) 

Here, Ai is the surface-type region A\ = {rG A s k/U(r, q) fl dAi, n 7^ for only one value of n face} 
where each r G A± is related to a unique face. Let us define Ai, n = 
{r G *4i/£/(r, <;) fl dAi <n 7^ for the n-th face} thus {.Ai,i, • • • , Ai }U , ■ ■ ■ } is a parti- 
tion of A\. Since r G .4i then r G .4i jn for some n face. Hence, the shape of U(r, q) fl .4 
is determined by the one dimensional position variable = d(r, dAi, n ) = z, and 

then, g(r) = gi(r) = g\ (r^) = g\ (z). In Eq. (130]) A2 is the edges-type region 
A2 = {r G A s k/U(r, q) fl <9.4i )m 7^ for two and only two values of m faces} where each 
r G A2 is related to a unique pair of faces. Let r G A2 and m 1 ,m2 be such that 
[/ (r, n dAi >mi 7^ and U (r, n ct4i im2 7^ 0, given that < 2q < £{A) the pair of faces 
1711,1712 converge to a common edge. Let us define A2, n = {r G A2I U(r,q) fl dAi, m 7^ 
for the two values of m faces that join in the n-th edge} thus {.4.2,1; • • • > ^2,n, • • • } is a 
partition of „4 2 - Since r G ^2 then r G ^4.2,n f° r only one edge (the n-th). Hence, the 
shape of U(r, q) fl A is determined by both, the dihedral angle and the two-dimensional 
vector that lies in a plane orthogonal to cM.2,n and goes from dA2. n to r, and thus, 
g(r) = g 2 (r) = g 2 (r^). A3 is the vertex-type region A3 = {r G Ak/^(r,?) H <9*4i, m 7^ 
for three or more values of m face} where each r is related to a unique set of faces. Let 
r E A3 and {mi, m 2 , m 3 , ...} be such that U(r,q) fl cL4i >mi 7^ 0, U(r,q) fl dA\^ m2 7^ 0, 
£/(r,^) fl &Ai im3 7^ given that < 2q < £(A) the set of faces {mi, m 2 , 7713, ...} 
converge to one vertex. Let us define A3, n = {r G ^.3/ £/(r, ^) fl dA\, m 7^ 
for the z values (with i > 3) of m faces that join in the n vertex} thus {.4.3,1, ■ • • , A3 >n , • ■ ■ } 
is a partition of A3. Since r G ^.3 then r G A3, n for some n vertex. Hence, the shape of 
U (r, q) fl .4 is determined by the three-dimensional vector that goes from dA3, n to r 
(and the angles that define the geometry of the vertex). Naturally, g(r) = (73 (r) = g 3 (r®). 
Care must be taken about dA m ^ n [the (m, n) element of OA] that should not be confused 
with d (A m ,n) [the boundary of the (m, n) element of the partition of .4] an irrelevant 
magnitude in PW. 
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Using the partition introduced for Ai, A2, A3 and the properties of g(r) we obtain 

I g{r)dv= [ 9i(r)dr+ ^ j g 2 {v)dv + ^ f g 3 (r)dr , (31) 

^•^ sk n faces ^ -^l.n n edges -^ 2 ' n n vertex ■'•^S.n 

where the right hand side term is a contribution to each vertex similar to that appearing in 
Eq. (|27p . Each of the edge's terms 

f g 2 (v)dr, (32) 

can be treated separately. Let us consider the family of orthogonal planes to the direction 
of the 72-th edge which intersects dA 2 ,n- Each of these planes cuts A 2 , n producing a slice or 
orthogonal cross section of A 2 , n - Let A 2n be the cross section defined by one of such planes 
that intersects A2, n in a region where both endpoints of the n-th edge are away to ensure 
that the shape of A^ n does not depend on any arbitrary choice. One can define £> 2 , n as 
the righted version of A2, n , which is the right (generalized) cylinder obtained by translate 
A^n along dA2, n , and extend the domain of g 2 {?) to B 2jn - The set B>2, n \ A 2 ,n contains two 
disconnected regions one around each endpoint of dA2, n (related to a given vertex). Each of 
these regions is (i3 2i „ \ A2, n ) n i with n' running over the vertex of A and (#2,71 \ ^2,n) n / — 
if edge n and vertex n! are not neighbors. Therefore 

/ g 2 (r)dr= / g 2 (r)dr- ^ / ^(r)dr , (33) 

J -M,n J &2,n n'vertex {&2,n\A2,n)„j 



with 



B 2 ,n JAi 



g 2 (v)dr = L n I g 2 (r®)dr® , (34) 



">2,n 

where is the two-dimensional coordinate of a point in A^ n and the last integral is a 
function of the dihedral angle at the n-th edge. Hence, for the central term in the right- 
hand side of Eq. (l3"Tj) we found 

J 92(r)dr= L n J ± g 2 (r^)dr^ - E E / g^)dr^ . 

n edges A2,n n edges -^2,n n'vertex n edge {&2,n\A2,n) n i 

(35) 

Now, to analyze each of the face terms in Eq. ( )3~Tj) 

</i(r)dr, (36) 

•^■1,71 

let us define B\ tn as the righted version of A\ lTl , which is the right (generalized) cylinder 
obtained by translate the face dAi, n along qz n being z n the inner (to A) normal unit vector, 
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note that Ai, n C B\, n . One can extend the domain of <7i(r) to all the semi-space that 
contains B\, n and includes dA\, n in its boundary plane to obtain 

f gi(v)dv=( 9l (r)dr- [ 9l (r)dr , (37) 

J Al n J 01 n J Bl n\A-\ r, 



with 



f 9l {r)dv = A n [ 9l (z)dz, (38) 

JBi, n JO 

where the last integral is independent of the involved face. Besides, £>i 5 „ \ Ai )H is a set of 
points distributed in the neighborhood of the relative boundary of A\, n i.e. rel<9 (^4i, n ) = 
Li n idA2,n' U„» dAs t n" with n' (n") running over the edges (vertex) that are neighbors of 
dAi )U . We make a partition of B\^ n \ A\^ n in regions corresponding to edges and vertex 

{• • • ' ^2,rm'> • • • ' ^3,nri"' • • } 

Bl,n\An= (J B V U 6 3,»». ( 39 ) 
n'edges n"vertex 

therefore 

I 9i(r)dr= J2 /, 9i(r)dr+ ^ / <7i(r)dr, (40) 

JB ltn \Ai, n n'edges J B 2,„„> n" vertex J B 'z, nn " 

where again n' (n") runs over the edges (vertex) that are neighbors of the n face. This 
partition is a simple extension of A^ n > and Aj,, n "i and it is warranted because the boundaries 
between A\, n A.i, n ' and >A.3,n" are built by pieces of planes, cylinders and spheres that can 
be trivially extended to such regions B\, n \ Ai jn £ A. For the n'-th edge we consider the 
family of orthogonal planes that cut B' 2nn , each plane producing a slice or orthogonal cross 
section. Let (B 1;n \ Ai^)^, be one of such cross section that cut cL4 2 ,n' in a region where 
both endpoints of the n' edge are away to ensure that the cross section does not depend on 
any arbitrary choice. For each n' edge we define B 2nn , as the right (generalized) cylinder 
obtained by moving (#i, n \ Ai >n )^i along dAi, n '- We note that B' 2nn , \ B' 2nn , contains two 
disjoint parts one around each endpoint of dA2, n ' which we assign to the corresponding 
vertex. Each of these contributions is (B 2nn , \ B 2nn ,) n „ w ^ n n ' running over the vertex and 
{&2,nn' \ ^2,nn') n // = if face n, edge n' and vertex n" are not all neighbors (taken in pairs). 
Therefore, for a given n face and n' edge we have 

f 9l (r)dr= [ 9l (r)dv- £ ( 9l (r)dv , (41) 

2,nn' 2,nn' 7l"vertex \ °2,nn' \ 2,nn f ) n tt 
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(42) 



with 

0i(r)dr = L n , / #i(z)dr (2) . 

One can verify that the edge contribution in Eq. (l42p for each of the two faces that meet at 
the n'-th edge is the same. By adding the contribution of all the faces, i.e. joining results 
from Eq. (ETJ) to Eq. (Ii2"| . one found 

[ 9i(r)dr = A P 9l (z)dz - 2 ^ L n , f ± gi(z)dr& 

ra faces ^- 4l ' n ^° n'edges J (Bi, m \Ai, m ) n , 

)dr® - [ 

JB' 



n"vert n face 



.n'edge l^2,nn'\^2,rm' J„« 



^i(2;)dr 



(3) 



3.nn ,/ 



43) 



where the m label (which is not an index) corresponds to any of both faces that meet at the 
n' edge. Putting all together, one obtain 



t = Vc + Ad + ^ L nC2,n + E ° 3 ' n ' 
n edges 



n vertex 



with the coefficients 



c 2 > 



Ai 



Co = G(r) for all r E A 

c i= f 9i(z)dz , 



(Bl,m\^l,m); 



g 2 (r®)dr® - 2 



g 1 {z)dr 



(2) 



(44) 

(45) 
(46) 

(47) 



c 3 



,» = / ^ 3 (r (3) )^ (3) - E / <fc(r (2) )dr( 3 > 

^•^3,n n'edges ^ ( S 2,n' Wn' ) n 



.n'edges ^ ( B 2,n"n' \ B 2,n"n' J n 



R' 

o,n" n 



(48) 



n"faces 

and g(r) = G(r) — c . ■ 
Theorem 2 (case d=2): For the case d = 2 Eq. (12"S"|) reads 

t = [ G(r)dv 

J A 

= Ac + Lc x + ^ c 2 , n , (49) 

n vertex 

where A is the surface area of A and L is the length of its perimeter OA. Besides, Co and 
Ci are constant coefficients which are independent of the size and shape of A, while c 2 ,n is a 
function of the angle in the ra-th vertex. 
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Proof (case d=2): For brevity, we only present such parts of the demonstration that 
differs from the case d = 3. Introduce a partition of A given by {A , -Ask} to obtain 

f G(r)dr = Ac + [ g{r)dr . (50) 

J A JA sk 

Consider the set Ask partitioned in terms of {Al,^}, 

Ak = AiUAz. (51) 

Ai is the line-type region Ai = {r G A s k/U(r, q) fl dAi >n ^ for only one n side} where 
each r G A\ is related to a unique side. Let us define A\ >n = {r G Ai/U(r, q) fl dAi )H ^ 
for the n-th side} thus {-4i,i> • • • , A\ jn , • • • } is a partition of A\. Since r G A\ then r G Ai jU 
for only one side (the n-th). Hence, the shape of U(r, q) (lA is determined by the one dimen- 
sional position variable r*- 1 -* = d(r, dAi, n ) — z , and then, g(r) = gi(r) = g\ (r^) = g\ (z). 
A2 is the vertex-type region A2 = {r G A s k/U(r, q) H dAi yin 7^ for two and only two 
values of m side} where each r G A2 is related to a unique pair of sides. Let r G A2 
and m 1 ,rri2 be such that U(r,q) fl &4i,mi 7^ and Z7(r, c) fl dA\ >m2 7^ 0, given that 
< 2q < 2(A) the pair of sides m^mj converge to a common vertex. Let us define i 2](1 = 
{r G A2I U (r, ij) fl dA\. m 7^ for the two values of m sides that join in the n-th vertex} 
thus {.A-2,1, • • • ,^2,11, • • • } is a partition of Ai- Since r G A then r G *42,n for only one 
vertex (the n-th). Hence, the shape of U(r, q) fl A is determined by both, the vertex angle 
and the vector that goes from dA<i,n to r, and thus, g(r) = g-z (r) = g2(r^). 
Using the partition introduced for ^4 1; .4 2 and the properties of g(r) one obtain 

f g(v)dr= [ 9i(*)dr+ ^ f g 2 (v)dr , (52) 

J A^ n gides tAl > n n vertex *'- 4 2 J n 

where the right hand side integral is part of the n-th vertex contribution to t, similar to that 
appearing in Eq. (I49D . To analyze each of the side terms in Eq. (|51|) 



f 9l (r)dr, (53) 

JAi,„ 

it is convenient to define B\ jn as the righted version of Ai >n , which is the rectangle obtained 
by translate the side dA^ n along qz n being z n the inner (to .4) normal unit vector, note 
that Ai, n C Bi tn . Let us consider the domain of #i(r) extended to all the semi-plane that 
includes dA\ >n and contains B\^ n . In terms of B\, n the Eq. (153]) can be written as 

/ gi(r)dr= g x (r)dr- gi(r)dr , (54) 

•J Al, n J Sl, n J Bl,n\Al,n 
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with 

/ gi(r)dr = L n / g\{z)dz , (55) 

where L n is the length of the n-th side and the right hand side integral is independent of 
the involved side. Furthermore, Bi >n \ Ai jU contains two disjoint parts one around each 
endpoint of dAi jTl which one assign to the corresponding vertex. Each of these contributions 
is (Bi >n \ Ai 7n ) , with n' running over the vertex and (Bi >n \ Ai iTl ) , = if side n and vertex 
n' are not neighbors. Adding the contribution of all the sides, joining results from Eq. (|54|) to 
Eq. (155]) and taking into account that each vertex contributes (with identical contribution) 
to a pair of sides it is found 

J2 I gi(r)dr = L T 9l (z)dz-2 ^ f gi(z)dr® (56) 

n sides*'" 41 .™ ^° n'vertex {Bi, m \Ai, m ) n i 

where the m label corresponds to any of both sides that meet at the n' vertex (and each 
side must be considered once). Putting all together, one obtain 

t = Ac + Lc x + ^ c 2,« > ( 57 ) 



n vertex 



with the coefficients 



c = G?(r) for all r G A , (58) 
ci = f gi(z)dz , (59) 



o 



C2,n 



f g 2 (r^)dr^-2[ gi(z)dr® . (60) 



Theorem 2 (case d=l): For the case d — 1 the Eq. (I2"6"j) reads 



G(r)dr 

Lc + 2ci , (61) 



A 



where L is the length of A while dA is composed by two separated points. Besides, cq and 
C\ are both constant coefficients which are independent of the size A. Note that, being A a 
straight line its shape is unique. Furthermore, in this case L = 2(A). 

Proof (case d=l): Again, we present such parts of the demonstration that differs from 
the cases d = 2,3. Introduce a partition of A given by {A Q , A s k} to obtain 
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f G(r)dr = Lc + f g(r)dr . (62) 
J A JA& 

with g(r) = G(r) — c and c = G(r) for all r G A . Now, ^4 sk — Ax — {^1,1,^.1,2} is a 
pure vertex-type region. Here, Ax,i and Axp are disjoint regions and each one corresponds 
to a given vertex or point, being dAx >n C Ax, n - Since both vertex have the same shape 
and the shape of U(r, q) H Ai >n is determined by the one dimensional position variable 
= d(r, dAx, n ) = z, then g(r) = <7i(r) = ^ (r^) = g x (z). Thus, one find 

f g{v)dv = 2 Cl = 2 r gi (z)dz. (63) 
JA Bh JO 

The Eqs. ( 1621 163|) can be re-arranged to obtain 

t = Lc + 2ci , (64) 

with 

c = G(r) for all rG A, (65) 
cx = f 9i(z)dz . (66) 







Finally we can broaden the field of application of Theorem 2 from P d to D P d . We 
define A G P^ if ^4 C ]R d , and ^4 is the connected union of count ably many closed convex 
polytopes, and 3 e G M>o such that VA G (0, e) and for every r G A the set int [[/ (r, A) R ^4] 
is connected. In this case it is sufficient to consider the extension of Eqs. (127|) and (|49p to 
the case of sums over countable many faces. Note that even when ^4 G P^ the theorem 
produce useful results only if £(A) is positive definite. 



A. Application of both theorems to confined fluids 

Now we turn the attention to a fluid confined by a polytope A G P (the fluid was described 
in Sec. [Til) an d analyze the functions rj : P — > M, : P — > R and : P — >■ R. In particular, 
to analyze r i; one must first consider the Theorem 1 and thus apply Theorem 2 to solve the 
integral in Eq. (1221) . We have obtained the following corollaries of the Theorems 1 and 2 
that apply to a confined fluid (relations are written for the case d = 3 but other values of d 
are easily included in our approach): 
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1. If a system of particles that interact via a pair potential of finite range £ is confined in 
a polytope A G P 3 such that £(A) > 2q with <^ taken from Eq. (1201) for some positive 
integer % > 1 then, the z-th cluster integral Tj : P 3 — > W is a linear function in the 
variables V, A and {L 1; L 2 , ...} (the length of the edges of ^4). It follows by replacing 
in Eq. (12 7p the magnitudes t, G with r,, E% and in Eqs. (jUJ) to (148 p the magnitudes 
t, G, Co, ci, C2 ir) , and C3, n with r i; z!&j, — i!a i; «!c|„ and i!cj„, respectively. Then, 
g(r) = -Ei(r) — Co (being Co = -Ei(r) with r G *4 ) and the expression of the i-th cluster 
integral is 

n l%\ = Vbi- Aa l+ J2 L nCln + E °ln > ( 6T ) 

nedges nvertex 

where the coefficients of b{ and <2j are independent of the shape of A while c^ n is a 
function of the dihedral angle in the ra-th edge and cj n is a function of the dihedral 
angles involved in the n-th vertex. Besides, all the coefficients depends on the pair 
interaction potential, the temperature and i. Finally, Eq. shows that b\ = 1 and 

a l = c l,n = C l,n = 0- 

2. If a system with iV particles that interact via a pair potential of finite range £ is 
confined in a polytope A G P 3 such that -2(^4) > 2^" with q taken from Eq. ( 1201) 
and i replaced by N, then its canonical partition function is polynomial on V, A and 
{Li,L 2 , ...}. It follows from the previous corollary, Eqs. ([2]) and (Q. From them the 
polynomials and Q can be explicitly obtained. 

3. If a system with at most M particles that interact via a pair potential of finite range £ 
is confined in a polytope A G P 3 such that £(-4) > ^ with ^ taken from Eq. (1201) and 
i replaced by M then, its (restricted) grand canonical partition function is polynomial 
on z, V, A and {L±,L2, ...}. This follows from the second corollary and Eqs. (JTJ)- (|3j) . 
From them the polynomial "Em can be explicitly obtained. 

For practical purposes one introduce the mean value for the edge and vertex coefficients 

C| = L- 1 Enedges L i C i,n aIld C I = ^ Envertex C i,n, respectively (with L e = Enedges L n and 

N v the quantity of vertex in A). Using this mean values the Eq. ( 1681) is re- written as 

nl%\ = Vb i -Aa i + L e ct + N v cJ . (68) 

with a\ = c\ = c\ = 0. 
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As a closing remark we can mention that for d = 1 it is easy to demonstrate that the 
conjectured relation (|21|) is true. Therefore, one found that n/i\ = biL — 2a« apply to any 
i < imax with i max = IntegerPart (L/q) + 1 [see Eq. fEE])]. 

B. Classification of g 

The three corollaries of Theorems 1 and 2 fix the structure of Tj, Zn and S, and provide 
explicit expressions that enables to obtain the coefficients of Tj. To make a further step in 
this direction we analyze the different scenarios for Ei(r) and g(r). At fixed r G A both 
functions depend on the set U (r, q) fl A throw the shape of U (r, q) fl dA. Once again, we 
mainly concentrate in the case d = 3. In terms of the partition for A introduced in Eqs. (128]) 
and ( 130]) r G ^4 =>- r G (for some fc = 0, 1, 2, 3) and then g(r) = <7fc(r). We have considered 
all the possible polyhedrons A with £(A) > 2q and observed that there exist countable many 
different cases for the geometry of U(r, q) D A and in each of this geometries g(r) must be a 
smooth function of r and the angles that describe the local shape of dA. In Tabled] we present 
such a classification for d = 3 and the case where at most three faces join in each vertex. 
There, the increasing complexity of the scenario with increasing k is apparent, there are 
only one case for k = 0,1, three cases for k = 2 and nine cases for k = 3. This behavior is a 
consequence of the number of independent parameters involved, which are 0, 1, 3, 5 (none; z; 

and (3; and and two -of the three- angles a) for k = 0,1, 2, 3, respectively. This table 
shows e.g. that for a given r< we must know three different expressions for g<i = g2(r ( - 2 \ a) if 
we expect to obtain the function c|(a) which enable to describe all the possible c\ n . Finally, 
Table |T] also apply for d = 1,2. For d = 2 we must restrict to < f < 2 and drop column 
v. On the other hand, the columns f and e should be renamed as: s (sides) and v (vertex), 
respectively. For d = 1 we must restrict to < f < 1, drop the columns e and v, and rename 
column f as v (vertex) . 

V. EXTENSIONS TO NON POLYTOPE SHAPE OF THE CONFINEMENT 

Once we have demonstrated the theorems and corollaries that reveal the shape of t(A) 
(and other magnitudes) for A G P with a characteristic length 2(A) > 2q, we are ready to 
relax some of the constraints on the geometry of A. In a first step we analyze the case of 
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Table I: Classification for the different shapes of U (r, q) f] A (dark region of the draws, blue online). 
The k index means that r G Ak and g(r) = guix^) with go(r) = Ei(r), the j index is introduced 
to enumerate the different scenarios for the shape of U (r, q) f] A at a given k. Columns f, e and 
v count the number of faces, edges and vertex, respectively, that intersects the ball centered at r 
with radii q. 



curved faces. Focusing on IR 3 , if one of the planar faces is replaced with a curved surface with 
constant and small curvature along it (cylindrical or spherical surfaces), what can we expect 
to obtain for the behavior of t(A)?. On the basis of the locality and translational invariance 
of Ei(r), the Eqs. ( 1221) and (1671) . and the ideas developed in the proof of Theorem 2 (case 
d=3), we expect that the coefficients corresponding to the curved face, their vertex and edges 
transform to functions of the principal radii of curvature of the curved face Rj, Rjj (and the 
involved dihedral angles). This functions of curvature radii could be expressed in powers of 
its curvatures Rj 1 and Rjj, being the independent term of each series that corresponds to 
the planar case given in Eq. (1671) . For several faces with small constant curvature one follow 
the same argument, in this case the edges (vertex) coefficients will depend on the curvatures 
of the faces that meet in the edge (vertex) and also on the involved dihedral angles. Related 
to this curvature dependence, in the case of a hard spheres (HS) fluid the cluster integrals 



r 2 and r 3 were obtained for a system confined by spherical walls 
evaluated for HS confined by cylindrical and planar walls 



10 



10Ml2|. while r 2 was also 



13]. In addition, the analysis 



of the HS fluid confined in the rf-dimensional space have provided explicit expressions of t 2 



in the cases of spherical and slit, confinement 



14L Il5| . Furthermore, in the case of square- 



25 



well particles T2 was evaluated for a spherical wall confinement [16]. In all these cases the 
obtained expressions coincide with the above discussed general picture. 

A second generalization of the theorem corresponds to the periodic boundary con- 
ditions which are frequently used in molecular dynamic and montecarlo simulations. 
We consider a set of three non-coplanar vectors {v 1 ,v 2 ,v 3 } that define the cell C = 
{r/r = X1V1 + X2V2 + X3V3 ,V0<Xj<l}, a parallelepiped, and the corresponding vectors 
of the Bravais lattice v = miVi + m2V2 + 1713V3 with m« running over all the integers. The 
translation of the cell over all the Bravais lattice vectors v tiles IR 3 . Note that systems 
with discrete translational symmetry in only two or one or zero directions (the last case 
corresponds to the absence of discrete translational symmetry) can be obtained in two ways; 
by fixing m 8 = in the definition of the Bravais lattice vectors for one or two or three of 
the integer numbers {mi, 7712, 7713} (which reduce the lattice dimension), or alternatively, by 
taking the large cell limit in one or two or three of the lengths {^1,^2,^3} with Vi = |vj|. 
We will take the second point of view to analyze the periodic polyhedron problem in the 
framework of the three dimensional Bravais lattice. Let us define the periodic polyhedron A 
by giving, a set of lattice vectors {vi, V2, V3} that define C and v, a polyhedron B C C and 
the periodic array of polyhedrons formed by all the copies of B over the lattice, Bi at . Fur- 
thermore, we assume that B\ a t belongs to a kind of generalized P space where the possibility 
of certain types of union of infinite convex sets (that consistent with a Bravais lattice) are 
allowed. We have strong evidence that shows that theorems 1 and 2 and their corollaries 
can be extended to include the case of a periodic polytope A. In this case the characteristic 
length £ (A) should be redefined. On one side, one should consider the characteristic length 
of the normal boundary of A, £' (A) = £(Bi at ), given by Eq. (1251) . that can be measured 
on a subset of B\ a t given by the cell and its 3 3 — 1 neighbors. On the other side, one should 
take into account the second characteristic length of A which concerns the geometry of the 
periodic cell, given by £" {A) = min (v 1, t>2, 1*2) • The minimum between both lengths £' (A) 
and £" (*4) is the characteristic length of the periodic polytope £(A). Due to at present, 
we have not a full demonstration of this generalization to periodic polytopes, it must be 
considered a conjecture. 

Other interesting extensions that should be analyzed in the future are: multi-component 
systems, system with non-spherical interaction potential, and external potentials that in- 
clude other non-hard wall-particle interactions. 
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VI. THERMODYNAMIC PROPERTIES 



In Sees. II III and IIVI it were derived new exact results on the statistical mechanics of con- 
fined fluids inhomogeneously distributed in the space without any particular spatial symme- 
try that applies both, to systems composed of few particles, as well as to systems composed 
by many particles. Recently, a thermodynamic approach that enables to study in a unified 



way both types of systems was formulated 



13 



17fl. We apply here this thermodynamic 



approach to study the properties of a fluid confined in a polytope in the framework of both, 
canonical and grand canonical ensembles. This approach rest in two basic assumptions 
made about the system: (1) it can be assigned to the system a free energy, this can be done 
through the usual free-energy -H- partition-function relation which is a well established link 
between statistical mechanics and thermodynamics, (2) under general ergodic conditions 
there exist an identity between: the mean temporal value of the properties of the system, 
those obtained from a thermodynamic approach, and those found from mean ensemble value 
(the ensemble that mimics the real constraints on the system must be considered). In the 
rest of this section we will find and discuss the thermodynamic relations, the equations of 
state (EOS), the relations between different EOS, and their meaning. 

For a system with a fixed number of particles, N, and temperature, T, one must adopt 
the canonical ensemble. In this framework, the connection between statistical mechanics 
and thermodynamics is 

F = -(3~ 1 \nQ, (69) 

where F is the Helmholtz free energy of the system. Besides, the thermodynamic fundamen- 
tal relations for F that follows from the first and second laws of thermodynamic through a 
Legendre transformation are 

F = U-TS , (70) 

dF = -SdT-dWr. (71) 

Here, S = —dF / 8T\n > j l is the entropy (the derivative is taken at fixed size and shape 
of the vessel), dW r = Py/dV is the total reversible work performed by the system on its 
environment and P\y is the pressure-for-work [13j|. For a general vessel A with a complex 
shape an arbitrary transformation may modify its size and shape. Thus, it is convenient to 
fix a given transformation for A and its parametrization. Using the A parameter one found 



27 



A — > A(X), dV — > dX dV/dX and the first EOS of the system 

dF 



Pw -~dx 



T,N 

It is interesting to note that once the expression for F(T, X) and Py/ are known an apparently 
thorough description of the thermodynamic magnitudes involved in Eqs. (JTOllTTj) is obtained. 
On the other hand, given that the system under study is an inhomogeneous one Py/ is not 
the pressure in the fluid, in fact, this simple observation shows that the obtained description 
of the system in terms of Pw alone is unsatisfactory because it said almost nothing about 
the properties of the fluid itself. 

To make further progress toward a more satisfactory thermodynamic description of the 
system properties, we assume that A G P 3 is a SB polyhedron and the confined fluid is 
such that the Corollaries [I] and [2] in Sec JIV Al apply. Thus, based on Eq. (|67|) we adopt 
the continuous measures X = {V, A, L, a} where bold symbols L and ot are short notations 
for the set of edges length and dihedral angles, a. Naturally, these measures can also be 
parametrized with A when it becomes convenient. Besides, from Corollary [2] and Eqs. ( 1701 
[71]) . we obtain the explicit form of F, S, U (from the statistical mechanical mean value 
recipe U = d/3F/dj3\N,A) an d Pw i n terms of X. We also found the following EOS 



dF 
dV 



dF 



dF 

7 ~ OA 

dF 



(73) 

T,N,X-A 



where the partial derivatives can be interpreted as the response of F to an small change 
of V (A, L n , a n ) while the other measures are kept constant. The Eqs. (1731 Ell) ma Y be 
seen simply as definitions for P, 7, T n and u n . In this context, each of these magnitudes 
is related to certain type of work, for example P is the work needed to change the volume 
of the system from V to V + dV (at constant A, L and a) and u) n is the work necessary to 
change the dihedral angle a n in da n (at all the measures X but a n constants, which affect 
both the edge and their vertex). However, the Eqs. (1731 \7M also suggest the meaning of 
P, 7 and Tn-, hence, P should be the pressure in the fluid, 7 its wall-fluid surface tension 
and T n the line tension corresponding to the n-th edge. In favor of this interpretation one 
observe that when the fluid develops a region with homogeneous density the pressure P from 
Eq. (173|) is the pressure of the fluid in this region. This fact can be demonstrated using 
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Theorems 1 and 2 and an approach similar to that developed in Ref.(l3|. Thermodynamic 
magnitudes P, 7, T n and u n depend on the adopted measures X but are independent of the 
adopted transformation. An equilibrium equation relates the different EOS, 

p-p w = q il+ ( qr T) + ( quU ) , (75) 

where the q coefficients are of geometrical nature q 1 = (^) , (qrT) = Enedges QT,n7n, 
M = E^ft-.nWn, with q T . n = (f)" 1 <fc and q^ n = (g The Eq. fl?5j) 
is very similar to the Laplace Equation for a drop that express the equilibrium between 
the inner pressure (in the liquid phase), the outer pressure (in the vapor phase) and the 
liquid-vapor surface tension. Hence, we say that Eq. Q75j) is a Laplace-like relation. For the 
case of g = where Eqs. fljj) and ([75]) have none sense one may define rfW 7 = — 7iycL4, 
^ w = ff It (df) 1 and ^ e Laplace-like relation gives 

lw ~ 7 = ( r rT) + (r w w) , (76) 

Where (r T T) = Enedges r T,nT n (V U UJ) = E n vertex^,n W «. with r T,n = (^f)" 1 ^ aIld 

r ^,n = (^x) 1 ■ Eq. (I76"|) makes explicit that the usual approach to the surface ten- 
sion which consist in measure the reversible work of doing a transformation at constant 
volume does not give the thermodynamic surface tension 7, even it measure 7^ for this 
adopted transformation. Now we focus on the case of some simpler shapes for A to obtain 
more explicit expressions of Eq. (175|) . For a convex polyhedron A which is also regular all 
its edges and vertex are equivalent (equal dihedral angles and edges length) , one can make a 
re-scaling transformation (at fixed angles) by multiplying the vectors that fix the position of 
each vertex by a parameter /c«l, In this case, for an initial state X = {V Q , A a , L D } we found 
V = 3k 2 V , A' = 2kA , L' n = L n , , g 7 = (q T T) = T (with T = T n for any n-edge) 
and (q^u) = 0. On the other hand, for the case of polyhedron A such that its complement 
A c is convex and regular (if A c is convex we say that A is corn-convex), starting with an 
initial state X = {V = — V a , A ,L }, we obtain V = — 3k 2 V q , A' = 2kA q , L' n = L n o , 
which produce a minus sign in g 7 = —^7- and {qrT) — ~ T ^l'y ■ By fixing k — 1 the Eq. 
(175|) for both, the convex and corn-convex regular cases, reduce to 

v--v*«-'w. +T k' (77) 

where sg gives the sign of dV d ^ (which is 1 or —1 for the convex and corn-convex cases, 
respectively). This equation applies not only to a cuboidal confinement, both for the fluid 
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in the cuboid and for the fluid surrounding it, but also, for the tetrahedron and other 
platonic solids (octahedron, dodecahedron and icosahedron) , evenmore it also applies to the 
cuboctahedron and the icosidodecahedron confinements (see [25] for geometrical details). It 
is interesting to note that for the above analyzed cases, where Eq. (177)) apply, Pw reduce to 
the mean pressure at the wall (wall-pressure) related with the mean density on the wall p w 
through the contact theorem 

Pw = /TVw • (78) 

Open systems with inasmuch M particles and T, fi -or z- and A fixed must be analyzed 
in the framework of the grand canonical ensemble. In this ensemble the connection between 
statistical mechanics and thermodynamics is 

fi = -/3~ 1 lnS, (79) 

while the thermodynamic fundamental relations that follow from the first and second laws 
are 

n = U -TS - fiN , (80) 
dQ = -SdT -dW T - Nd/i. (81) 
By introducing a parametrization for the cavity transformation one find 



w 



on 



dV 
~d\ 



-i 



(82) 



For the case of A G P 3 a SB polyhedron and a fluid such that the Corollaries [T] to |3] apply, 
we use the same measures X adopted above. Besides, from Corollary [3] and using Eqs. (Tj[|) 
to (J7J) we obtain the explicit form of Q and Pw in terms of X. The expressions for S, U and 
the EOS in this GCE are easily found by replacing F with Q in Eqs. (172)) to (171)) . For the 
other EOS we obtain 



P 



dQ 
dV 



7 



T,z,X-V 



dQ 



dL n 



on 

dA 
dQ 



U) r , 



T,z,X-L n 



T,z,X-A 



T,z,X-a n 



(83) 



(84) 



As can be seen, several results are similar to that found for the canonical ensemble case and 
thus we only draw about some features of this GCE study. For example, from the EOS one 
obtain the equilibrium relation 



(85) 
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this Laplace-like relation is the same that Eq. (]75j) and again P is the pressure in the fluid, 
which coincides with the pressure in a region with homogeneous density [see details in Eq. 
(175|) or Eqs. ((771 ITS]) for the specific case of regular polygons]. Other usual thermodynamic 
and statistical mechanical relations that are valid even for the case of inhomogeneous systems 
with a finite M value are 

d(3tt 



N = (N) = -z 



dz 



T,A 

dN 



(J, 



(86) 



(87) 



where a N , the standard deviation in the mean number of particles N, quantifies its sponta- 
neous fluctuation. Note that derivatives with respect to z are taken with the fluid in a fixed 
region A (all the X measures fixed). 



A. The low z expansion 

To study the thermodynamic behavior of the confined open fluid in the low z regime 
one express Eq. (179]) as a power series in z [see also Eq. fl3])]. We introduce the vector of 
coefficients bj = (b{, dj, c| 1; ■ ■ • , cJi, • • ■ ) which depends on ex., and the vector of extensive- 
like magnitudes Y = {V, A, L, l v } with l v the vector of iV v components all of them equal to 
one being iV v the number of vertex. Hence 

Tj r = Y-h l , (88) 

an expression equivalent to Eqs. (1671) and ( 168]) . Using Eq. ([9]) we can write 

M 

W = ~J2^ zi = ~I> iy -bi + OM+tiz) , (89) 

i>l ' i=l 

which is exact to order M. Notably, Eq. (189]) shows that to this order Q is linear in the 
extensive-like magnitudes in Y, on the other hand, Eq. (168]) shows that to this order Q is 
also linear in V, A, L e and iV v . To the same order we found 

M 

N = J2iz l Y U, (90) 

i=i 

AI 

c = E jVF ' b - ( 91 ) 

i=l 
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Therefore, in Eqs. (l89|) -( !9T|) and to order M one can separate each component. For example 
the Eq. p| is 

tt = tt b + tt s + fi, + tt p , (92) 

with the volumetric part of the grand-potential and the grand potential of surface, 
line and points Q s , Qi and Q p , respectively. Each of they correspond to the sum of terms 
proportional to V, A, L e , N v ; and is trivially related with a grand potential density (per 
unit volume, area, etc.) 

M 

= /3" 1 ^^, (93) 



p _£h on 
b ~ v " dv 



■;T,X-V 



i=l 



7o 





dQ 
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~ dA 
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' oo 
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" L e " 






V - 


= iVv =; 



M 



z,T,X-A 



Z dj 



i=2 



-1 



M 

'E 

i=2 
M 

E 

1=2 



z l cj 



(94) 
(95) 
(96) 



Note that the right hand side equality in Eqs. ( 1931) to ( l96i) is up to order M in z. Here, Eqs. 
fl93l) and ( 1941 coincide with work terms definitions ( IHBl while the Eqs. (1951) and fl96l) may be 
understood as mean works. From the right hand side term of Eq. (193]) one find that P b is the 
pressure of the bulk system (with the same T and zY. This series was studied by Mayer and 
others (see Ref. [18] and a more complete survey in |6, p. 122]) in their approach to the virial 
equation of state for the bulk system. On the other hand, from the right hand side term 
in Eq. ( 194"1) one find that is the (wall-fluid) surface tension for the bulk fluid in contact 
with an infinite planar wall. This series is consistent with that obtained by Sokolowski and 



Stecki 19|, a question that will be briefly discussed at the end of this section. Furthermore, 
Pb and clearly corresponds to P and 7 from Eq. (183]) but also to the CE Eqs. (173]) . 



which validates our thermodynamical approach. The additional work terms, to order M in 
z are [from Eq. ( 184]) ] 

M 

Tn ^ ^ Z C i,n 1 



1=2 



M 



1=2 



<9c? 



+ 



dc, 



t.p 



+ 



dc. 



(97) 
(98) 



z,T,X-a n 
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where p and q are the vertex at the endpoints of the n-th edge. Eq. (l97j) corresponds to 
the line tension on the n-th edge. Thus, one can prove [using Eq. (168]) ] that To = (Tn) = 
^e 1 Snedges TnL n , i- e -> Too is a mean-work term (the work needed to increase the edges length 
in dL e of an edge with mean properties) it is the mean line tension. On the other hand, 

<9(iV» 

-i- 

z,T,X-a n 



T dT n 



+ da n 



(99) 

z,T,X-a n 



In terms of the densities the Eq. ( 192]) transforms to 

Q = -VP b + A loo + L e Too + N v u , (100) 

which is a demonstration that the grand potential is a homogeneous function of the so 
called extensive variables V, A, L e and N v (at least to order M in powers of z). Note 
that this is a central assumption in the thermodynamic and statistical mechanics theo- 
ries of macroscopic systems composed by many particles and spatially distributed follow- 
ing strong symmetries. On the opposite, here it is demonstrated for an inhomogeneous 
system, may be composed by few particles, without particular translational or rotational 
spatial symmetries. Naturally, explicit expression for several excess grand potential can 
also be obtained, for example the (over-bulk) surface excess grand-potential density is 
7 = (fi + VP b ) I A = (3' 1 Y;T=2 ^ (a, - c\L e jA - cJN v /A) which is a rough estimate of the 
surface tension. Introducing the A parametrization for a general vessel transformation we 
found the Laplace-like equilibrium relation 

P b -P w = g 7 7oo + (q T T) + (qM , (101) 

equivalent to Eq. (185]) . but in this case we found the expression and meaning of each term 
to order M in z. Again, for the particular cases analyzed below Eq. ( 176]) we obtain 

2A L 

(P b - P w ) x sg = 7oo — + ^ ^ , (102) 

where Too = T n = T is a function of the dihedral angle [details are given below Eq. (177]) ]. 
This Eq. suggest the following procedure to obtain information about 7^ and Too- Consider 
7 defined below 

7 = ^g x (A - P w ) = 7oo + Too > (103) 

all the magnitudes between the equal signs are simple to measure in a molecular dynamic 
simulation of the open system. In particular Eq. (178]) . which remains valid in the GCE, 
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provides a simple way to evaluate Pw- Hence, one can fix the vessel shape, T and a small z 
value, and then do measures of 7 along simulations for several different sizes of the vessel. 
Plotting 7 vs. |j and fitting the points with a linear regression we obtain from the ordinate 
and the slope 700 and 7^., respectively. If we repeat the procedure for different values of z 
we obtain a table of values for 700(2;) and Too(z) which enable to evaluate the coefficients 
and c| through a second fit in this case with a polynomial function. 

The splitting of Q in terms of its components to order M in z, given in Eqs. fl92|) and 
(llOOp . can also be done for Eq. (190]) . It gives 

M 

{p - Pb ) V = AT = AT + L e T e + N V T V = iz* (-M + L ec e i + iV v cJ) , (104) 

8=2 

where p = N/V is the mean number density, pb = Y^i=i iz z bi is bulk density, T is the 
total effective adsorption, while the area, edge and vertex adsorptions are T, T e and r v , 
respectively. On the other hand, the same procedure applied to the fluctuation in Eq. (19~TT) 
produce 

M 

a N 2 - o£ = As = As + L e s e + N v s v = ^ i 2 z i (- ai A + L & c\ + N v a[) , (105) 

i=2 

with a 2 = V Yliti i 2 z % bi the fluctuation in the number of particles in the bulk, s an excess 
effective fluctuation per unit area while s, s e , and s v are the components of surface-area, 
edges length and vertex of the number density fluctuation. Particular attention deserve such 
polytopes for which £(A) is infinite. They are the unbounded polytopes with one face, or 
one edge, or one vertex. In this cases M can take any positive integer value and thus it be 
made as larger as one wishes. Therefore all the series in powers of z given between Eq. ( 1891) 
and Eq. (11051) are valid to any order. Besides, the linear decomposition of several functions 
in its basic extensive measures Y (A) is also exact to any order. In particular, the grand 
potential becomes an homogeneous function of the extensive measures [see Eq. (llOOp ], 

The structure of Eqs. (I92H961 11001 H04p and fl 1 5 [) show that the relations between the 
functions 

(f3P b - z, p b - z, cxb - z) and 6j>i, (106) 

(where (3Pb — z, Pb~z and o 2 — z are excess magnitudes with respect to the ideal gas system) 
is the same that the relation between the functions 

(/?7, -f , -s) and di - c^LeA' 1 - cJNyA' 1 , (107) 
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while the same apply to these other sets of functions and coefficients 

G^Too, -r, -s) and a i; (108) 

(-pTao, r e , s e ) and c\, (109) 

(-/3z/, T v , s y ) and cj. (110) 

These simple observations show for example that — /37oo(r) and the dependence of (3Pb — z 
on pb — z is the same except that — a« in the first case must be replaced by hi in the second 
case. 

Trivial series manipulation (inversion and composition) enables to obtain explicit expres- 
sions for the power series (to order M) of P b (p), Pb(pb), loo(Pb), l(p), 7( r ), 7( s e), F(Pb), 
Too{p)i %o{Pb) and others. Several examples of these series discussed in terms of the formal 
symmetries given in Eqs. fllQ6p to (IllOp are explicitly evaluated in Appendix |A] For fixed 
M we have verified that the series for Pb{p b ) (to order M) is the usual virial series (see for 
example Ref. Jfl), while 7oo(p&) and T(pb) reproduce the approach of Bellemans j^Pj and 
subsequent improvements of Sokolowski and Stecki [19] which have developed a complete 
theory of virial expansion for surface thermodynamic properties for the case of a fluid in 
contact with a planar infinite wall. As far as we known the series representation of any other 
of the discussed functional relations were never analyzed before. 



VII. FINAL REMARKS 

The obtained polynomial structure of the CE and GCE partition functions of few- and 
many-body inhomogeneous fluid-like systems confined by polytopes is interesting in several 
ways. It shows that the statistical mechanical properties of many of these fluid-like systems 
with a fixed number of particles can be exactly described if the involved coefficients (which 
depend on T) are known. This consideration also applies to open systems if we conveniently 
restrict the maximum number of particles considered. The key point for these results was the 
analysis of the geometrical properties of the reducible cluster integrals for inhomogeneous 
systems. It is interesting to mention that from the decade of 1950 until present this integrals 



(introduced by Mayer and Mayer 2l| for the study of homogeneous fluids) were not the focus 
of much interest and were almost ignored in the study of inhomogeneous fluids. On the other 
hand, given that our approach to Theorems 1, 2 and their corollaries is constructive we found 
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r 1 (T 1 ) 

Figure 2: The size vs. number domain for the system of ./V particles in a cube with side L. Below 
the dotted curve the Corollaries 2 and 3 do not apply and thus the exact structure of the partition 
functions is unknown. 

explicit integral expressions for the coefficients. In a work in progress we evaluate some of 
the unknown lower order coefficients for a system of hard sphere particles. 

The results for the partition functions were complemented with a statistical mechanical 
and thermodynamical recipe that, based on the linear decomposition of cluster integrals 
enable to find the exact properties of the confined inhomogeneous fluid (for both the open 
and closed systems). This fact is notable because the studied confined inhomogeneous 
systems, which are far away from the many particles in large volumes condition, have not 
any particular spatial symmetry and thus the very basic assumptions of the usual approach to 
thermodynamics theory are not satisfied. Furthermore, it was shown that in the low density 
(and low z) regime the obtained grand potential is a homogeneous function of the extensive 
variables, which is a consequence of the linear decomposition of the cluster integrals in their 
extensive components (for d = 3 they are volume, boundary surface area, edges length and 
vertex number). Even more, it was found a generalized version of virial series EOS, that 
provides expressions for pressures, wall-fluid surface tension, line tension, adsorption, etc. 
All this findings strongly suggest that it is possible to build a wider scope formulation of 
both, statistical mechanics and thermodynamics, theories. These generalized formulations 
should concern to systems with or without symmetries, involving from few- to many-bodies, 
and constrained to regions of any size. 

In Fig. [2] it is represented the domain where the CE (GCE) partition function is a 
polynomial for the case of iV particles (at most iV particles) confined in a cube with side 
L. Broad straight lines show this domain for Zj^ and S^r for different (fixed) number of 
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particles. Each line extends from the infinite dilution limit L^ 1 — > to the maximum 
density for which the Corollaries given in Sec. IIVAI apply. This maximum density condition 
is given by AT -1 = (in accord with Eq. (12T|) and Sec. IIV Ap and is drawn in dotted 

line. The large systems limit corresponds to V — > oo, A" — > oo and coincides with the origin. 
However, to analyze the thermodynamic limit one must reach large systems by following a 
constant density path, which correspond to a dependence A^ _1 oc L~ 3 . In the figure it is 
plotted for a small value p = 0.2 £ -3 with a dashed line. It is clear that our results about CE 
are insufficient to analyze the thermodynamic limit because for large volumes (small abscissa 
values) the highest density under which the structure of the partition function is known gives 
a linear relation AT -1 ~ L -1 . On the basis of this result is the nature of Theorems 1 and 
2 which show that the end of validity of Eq. (I6T|) relates with the existence of a cluster 
configuration that is capable to percolate the cavity A. 

An attempt to analyze the polytopes A such that 2<j > £(A) shows two folds. On 
one hand, that for some of these polytopes the Eq. ( 1671) remains valid (as well as, apply 
the conclusions of Corollaries 2 and 3). In this case, appears a less restrictive condition 
2c' < £(A) with c' < c that replace 2c < £(A). On the other hand, that the end of 
applicability of Eq. (167|) relates with the non-universal behavior of due to the existence 
of a non-analytic term [which is identically null for 2c < £(A)\ that depends on the shape of 
A in a more complex way. The study of this term for some simple confinements is interesting 
because it may enlighten how to describe the system properties in the thermodynamic limit. 

The core of the formal result presented in this work relies on Theorems 1 and 2. As was 
mentioned, Theorem 2 for d = 1, 2, 3 and the conjectured generalization to M. d given in Eq. 
( I26|) strongly resemble well known results of integral geometry. In particular, the expressions 
given in Eqs. (126| [27] H9|) and f l6Tj) look like the combination of Hadwiger's characterization 
theorem and The general kinematic formula (Theorems 9.1.1 and 10.3.1 in pp. 118 and 153 

n 

of [9], respectively) as it were applied to a polytope and a ball. Even that, several differences 
can be underlined. On one hand, in PW it was not established if the functions G(A, r) and 
Ei(A, r) are valuations (also called additive functions) or not. On the other hand, the 
integration domain of the integral solved in PW is not the complete space. Furthermore, 
the question about the continuity or monotonic behavior of both functions on P (see p. 153 
in Ref. Q and pp. 211 and 253 in Ref. [7]) is open. Despite these differences, the connection 
of our results with various formulae from integral geometry and convex bodies, theories is 
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evident. For example, some parts of the presented formulation resemble to Steiner's formula 
for the measure of the Minkowski sum of a convex polytope and a ball (see Theorem 9.2.3 

n 

in p. 122 of [9]), while the solved integrals are similar to certain integrals over functions that 
depend on the distance between a point and the boundary of a convex set. (p. 132 of |9j and 
p. 258 of [7]). These connections should deserve a deeper analysis. 
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Appendix A: Series Expansions of the EOS 



We present here the expressions for several series expansions of thermodynamic magni- 
tudes. In the following Eqs. the series were truncated to the next order that were written. 
In the derivations we make intensive use of relations expressed in Eqs (I106p to (IllOp . From 
Eqs. fl94]) and (TT04j) we found 7^ (p b ) 

filociPb) = a>2pl + {a 3 - 4a 2 b 2 )pl + (a 4 - 6a 3 6 2 + 20a 2 bj - 6a 2 6 3 ) pi . (Al) 

Besides, —(3Too{pb) is identical to the right hand side of Eq. (jAl|) but with the replacement 
di — > c\. The same applies to — fiv{pb) with the replacement dj — > cj and also, to [i^f{p b ) 
with di — > ai — c^L e A~ 1 — c' i N w A~ l . From Eqs. (1941) and (I104p we found the relation between 
surface tension and adsorption 

„ , v 1 a 3 F 3/ 2 (-9«i + 8a 2 a 4 ) T 2 3 (63a| - 96a 2 a 3 a 4 + 32a^a 5 ) T 5/2 

Ploo{ ] ~ ~2 + 4v ^(-a 2 )3/2 + 32^1 + 256 v / 2(-a 2 )9/ 2 ' 

(A2) 

which shows a non-analytic dependence between fluid-wall surface tension and the surface 
adsorption in F = 0. Besides, — P%o(T e ) is identical to the right hand side of Eq. (1A2D 
but with the replacement dj — > cf, while the same applies to — (3v(T v ) with the replacement 
di — > cj and also to /37(F) with a« — )■ <2j — c^L e A _1 — cJA v A _1 . Furthermore, from Eqs (1941) 
and (I105P we found 

5 3 a 3 3/2 1 (81aj - 64a 2 q 4 ) g 2 

The same general ideas applied above in relation with Eqs. (LAI I) and (lA2j) enable to obtain 
e.g. —l3Too{.s e ) and —j5v{s y ). Fluctuation and density in the bulk, taken from Eqs. ( 11051) 
and ( I104p . relates by 

a 2 b = Vp b [1 + 2b 2Pb + (-8b 2 2 + 663) pl + {m 3 2 - A8b 2 b 3 + 126 4 ) p 3 b ] , (A4) 

while to obtain a 2 N we must replace p b — >■ p and Oj — > 6, — aiAV~ l + c^L e \^ _1 + cJiV v V -1 , in 
this case V p = N. On the other hand, s(p b ) is 

s(p b ) = Aa 2 p 2 b + (9a 3 - 16a 2 6 2 )p^ + (l6a 4 - 54a 3 6 2 + 80a 2 6 2 - 24a 2 6 3 ) p 4 b , (A5) 

while, s e (p b ) can be obtained through the replacement Oj — > cf, while s v (p b ) can be obtained 
through the replacement ai — > cj. To obtain s(p), we replace p b — >■ p and a« — )■ a» — c\L e A~ x — 
cJN w A- 1 in Eq. flA5l) . 
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